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Abstract 
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We find an exact analytic solution for the time evolution of a three Dirac 
neutrino system adiabatically oscillating in matter, constructing explicitly the 
relevant 3x3 mixing matrix in matter. Using this result we investigate the 
solar neutrino data in a scenario where the neutrino masses are such that 
mi < 771,2 ^ iTT'S, taking into account several phenomenological constraints on 
neutrino mixing angles and masses. A solution of the solar neutrino problem 
for large values of the parameter dm? = — m\ which are not usually 
associated with a resonance is found. This is an essentially three-generation 
effect. 

PACS numbers: 14.60.Pq, 96.60.Kx 
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Typeset using REVT^ 



Most of the analyses of the neutrino oscillation hypothesis assume that this phenomenon 
involves only two generations of neutrinos. It is difficult, however, to justify why oscillations 
would not involve also the third family. From the conceptual point of view, mixing and 
mass parameters required for three generation oscillations are not different from mixing and 
mass parameters that appear in this same phenomenon involving only two generations of 
neutrinos. Why would such parameters related with the third family vanish and not those 
ones leading to two neutrino oscillations? In our opinion, the two generation analysis is just 
an indicative approach to the more realistic three generation case. 

The larger number of mixing and mass parameters in three neutrino oscillations can be 
quoted as a difficulty to approach this scenario. There are three mixing angles, one phase 
and two independent mass parameters that are, in principle, free parameters. Due to this 
fact, when three generations are considered in the literature [||], some assumptions have been 
made to restrict this parameter space resulting that either only two family transitions are 
effective or the parameters are fixed arbitrarily. 

In this letter we find an exact analytic solution for the time evolution of a three Dirac 
neutrino system adiabatically oscillating in matter, constructing the three dimensional mix- 
ing matrix in matter, to calculate the electron neutrino survival probability and compare it 
to the solar neutrino data. 

Assuming the minimal extension of the standard electroweak model when only three 
right-handed neutrino singlets are introduced to generate Dirac neutrino masses, a mixing 
can then occur among the three lepton fiavors. Therefore, neutrinos produced in weak 
processes are in general linear combinations of the mass eigenstates: = J2i Vai^^i (ct = 
e, /i, r; i = 1, 2, 3), where 



ceCfB 



V 



CgCj — SgS^S/s 



S^Cp 



(1) 




— CqS^ — SgC^Sp 




We have set to zero the CP violating phase in Eq. (1). 
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The matter effects for the generahzed case of three generations are described by the 
time evolution equation H] 
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where 
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(2b) 



and A = \/2GFNe{t), with Gp the Fermi constant, Ne(t) the electron number density in the 
region reached by the neutrino at the instant t. 

In order to write the neutrino survival probabilities when matter effects are present, we 
would like to have the exact solution of the three coupled differential equations given by Eq. 
(0). This has proven to be very difficult to obtain and only approximate solutions has 
been achieved up to now. In the following we will obtain the explicit form of the relevant 
3x3 mixing matrix in matter. After that we will be able to construct an exact solution 
for the three neutrino time evolution equations assuming that the neutrino propagation is 
adiabatic everywhere. In this case the problem is reduced to diagonalize the matrix H in 
Eq. (^. We obtain the characteristic polynomial of H 



+ 3aX'^ + 36A + c = 0, 



(3) 



where 



3a = -Tr H, 3b = + H"^^ + H^^, c= -det H 



(4) 



and denotes the minor of the Haa elements. We will not write explicitly the elements 
Haa', «, a' = e, /i, r since they can easily be obtained from Eqs. (|I]) and (^. The eigenvalues 
are 



'0 + 2n7r' 

A„ = 2^/—h cos ( ) — a, n = 0, 1, 2; 



h = b — a^, g = c — Sab + 2a^, cos B 



The respective eigenvectors are 



where 



Vai = ^, a = e,/i,r; z = 1,2,3, 

Oi 



5,1 = (ifee - \o){H^^ - Ai) - Hl^, 5i = \{6elf + + {6rlf 



(5a) 

(5b) 



(6) 



(7a) 



5r2 = {Hee - A2)5^1 - -^e/x^el, ^2 = (5e2)^ + {Sfi2 f + (5r2) 



(7b) 



(7c) 



Je3 — ^Tl^i^a 

Therefore the phenomenofogical eigenstates can be written in terms of the matter eigen- 
states Ua = J2iVaii^i, {i = 1,2,3), where the matrix V can be read from Eqs. (H) and (^ 
and parametrize in terms of the mixing angles in matter as the matrix in Eq. (|l|) but with 
— s> 0, 7 — * 7 and (3 ^ (3. It is trivial now to write down the averaged adiabatic survival 
probability of finding a z/^ produced at the point Xq inside the sun and detected at the point 
X in the Earth's surface 



Pu^^u^=Y.\yo.^{Xo)?\V.^{x)\^ 



(8) 



The solar neutrino problem has been confirmed by many experiments. In the following 
we will consider experimental data from Homestake(H), Kamiokande(K) and Gallex(G) 
experiments [QQ. 

Neutrinos produced in different reactions have different energies. While ^Be neutrinos 
are almost monochromatic [§, neutrinos produced in other source- reactions have different 



energy spectra |^ which have to be considered since, as we will see in the following, the 
survival probability of the solar neutrinos is sensitive to their energy E or their momentum 
p. It is necessary also to take into account the value of the solar matter density at the 
neutrino creation point xq. We use the solar matter distribution calculated through the 
Standard Solar Model which is tabled in Ref. |^. Notice that in Eq. (^, since neutrinos are 
detected at the Earth's surface, the matrix elements at the point x are essentially those of 
the vacuum mixing matrix. 

We can compare the theoretical neutrino flux {(pth) calculated from the Standard Solar 
Model 0] with the observed flux {(pexp) measured by each experiment. The ratios R = 
(pexp/^th are given by i?(H) = 0.28 ± 0.04, R{K) = 0.49 ± 0.12,and R{G) = 0.66 ± 0.12 [|. 

Considering the neutrino oscillations, these ratios can be calculated for each experiment. 
We take into account only the main source reactions of solar neutrinos which are sensible to 
each specific experiment. For Homestake we have 

i?(H) = 0.78P"{^B) + O.UP^CBe) + OMP^C^O). (9a) 

The neutrino flux measured by Kamiokande facilities is not merely the electron neutrino one 
since detector electrons will interact with other neutrino flavors via neutral currents. For 
energies involved in the solar neutrino experiments, the z/e-electron scattering cross section 
is about seven times larger than other neutrino flavor (z/^-electron and z/T--electron) cross 
sections. Hence, for Kamiokande, taking into account these neutral current effects we have 

i?(K) = P^(8E) + ^[1-P^^f5)]. (9b) 

Finally, for Gallex 

R{G) = 0.2QP^{^B) + O.llP^CBe) + 0.05P^C^O) + OMP^{pp). (9c) 
In Eqs. (H) we use the notation 

P-\X)= f''{E)P,^^,^{E,Sm\e,x,Xo). (10) 
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J = H,K and G for Homestake, Kamiokande and Gallex; X denotes the particular source 
reaction of solar neutrinos and Pu^^ua is given in Eq. (^. The threshold energy for each 
one of these experiments and the energy spectrum of neutrinos produced in reaction X are 
denoted by Eff^,^^ and f^{E), respectively. The spectra function f^{E) are given in Ref. [§. 

As we said before, in Eq. (j^) there are still too many free parameters: three vacuum 
mixing angles and three masses. Hence, it is necessary to take into account the constraints 
from other physical processes to fix some of the neutrino parameters before taking into 
account the solar neutrino data. The neutrino masses and mixing angles for the case of 
TTLi ^ 1712 ^ ^3 have been determined in Ref. using r leptonic decays, pion decays, Z° 
invisible width and end-point data from r decay into five pions and assuming world average 
data for the ratio Gr/G^. Assuming the above mass hierarchy, the lower masses mi, m2 and 
one angle 6 remain undetermined, but ~ 165 MeV, 11.54° < f3 < 12.82° and 7 < 4.05°. 
Thus, we have one mixing angle 6 and two lightest neutrino mass difference 6m'^ = m\ — m\ 
to be determined in neutrino oscillation processes. 

Using Eqs. (|^) we have investigated the compatibility regions of the three solar neutrino 
experiments. We have considered the region 10"^ eV^ < (5m2 < lO^eV^ For 10"^ - lO^eV^ 
we have solution for either Homestake and Kamiokande or Gallex and Kamiokande but not 
for all of them at the same time. The only two regions fitting the three experiments are 
shown in Fig. 1 at 95 % C.L., we have restricted the range of 9 to be bellow tt as the survival 
probability is symmetric in tt — They correspond to: 3.5 x 10^ eV^ < dm? < 5 x 10'' eV^ 
with 0.9 < < 1.2 radians and 2 x lO'^eV^ < ^m^ < 2 x 10^ eV^ with 1.4 < ^ < 2.2 radians. 
In fact we have found small regions of allowed values even at 68 % C.L.. 

Some remarks are in order. With the parametrization of the mixing matrix used in this 
work, the resulting survival probability of the electron neutrino is not sensible to the angle 
7. Setting /5, ma, and A equal to zero, we recover the vacuum solution in two generation 



case On the other hand keeping A 7^ we recover the MSW solution in two genera- 
tions. In this case we have solution for 5m?' ^ 10~^eV^ and sin^ 2^ ^ 1, but only for two 
experiments at a time, what is consistent with the fact that we are considering only the 



adiabatic solutions [|TT] . 



General analytical descriptions of three generation neutrino oscillations are far from 
transparent. We have found a solution to the solar neutrino problem in a region of the rele- 
vant parameter space which is not usually associated with resonances in the three neutrino 
evolution equations This is a completely new feature compared with the two neutrino 
MSW solution where only in the vicinity of a resonance we can expect to find a dependence 
of the flavor survival probability on the neutrino momenta. In that case, if we do not have 
a resonance for the values of the oscillating parameters, all neutrinos undergo the same 
survival probability and it is impossible to conciliate all solar neutrino data. 

Here we have looked for such dependence of the flavor survival probability on the neutrino 
momenta investigating the behavior of the relevant mixing angles in matter considering the 
momentum range of neutrinos produced in the sun. While pp-neutrinos have momenta not 
larger than 0.44 MeV, ^5-neutrinos can present larger values of momenta up to 15 MeV. 
In Fig. 2 we show the results of this analysis. We observed that for values of the mass 
difference Sm^ and the vacuum mixing angle 6 found to be relevant for the compatibility of 
all solar neutrino data (see Fig. 1), a strong dependence of the mixing angle in matter 6 on 
the neutrino momentum p is present. 

Summarizing. We have found an exact analytical solution for the adiabatic transition 
probability in the three generation neutrino case. Applying our solution to an example 
where one of the masses is rather large, of the order of several MeV, we have shown that 
the mixing angles in matter strongly depend on the neutrino momentum when we consider 
the range of momenta physically interesting for the solar neutrinos. Therefore a solution to 
the solar neutrino problem can be achieved. 
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FIGURES 

FIG. 1. The compatibility region of the three solar neutrino experiments in the 5w?-9 plane at 
95% C.L. . 

FIG. 2. sin 2^ (^20) ^ ^ function of neutrino momentum p for dm? = 1 x 10^ eV^ and several 
values of 0. 
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